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Abstract 


In this study the elasticity problem for an infinite orthotropic flat plate 
containing a series of through and part-through cracks and subjected to bending 
and tension loads is considered. The problem is formulated by using Reissner’s 
plate bending theory and considering three dimensional material orthotropy. The 
Line-spring model developed by Rice and Levy is used to formulate the surface 
crack problem in which a total of nine material constants have been used. The 
main purpose of this study is to determine the effect of material orthotropy on 
the stress intensity factors, to investigate the interaction between two 
asymmetrically arranged collinear cracks, and to provide extensive numerical 
results regarding the stress intensity factors. The problem is reduced to a 
system of singular integral equations which is solved by using the Gauss- 
Chebyshev quadrature formulas. The calculated results show that the material 
orthotropy does have a significant effect on the stress intensity factor . 



Chapter 1 
Introduction 

In recent years the increasing use of composite plates in many engineering 
structures and especially in the aero-space industry has brought up the need for 
more intensive stress analysis of anisotropic materials. Up to now no completely 
analytical solution of the related failure problems has been available due to the 
inherent difficulties involving the stress analysis and material characterization. In 
many cases the failure is attributed to the growth of cracks or crack-like flaws 
which exist in the structure. During the past two decades many investigators 
have studied the stress state in the immediate neighborhood of the crack tip 
since the local fracture of a structure appears to be governed mainly by this 
stress field. The singular behavior of the stress state near the crack-tip as 
characterized by the stress intensity factor depends on the magnitude of the 
external loads, the configuration of the body including the crack size and shape, 
and the material properties. Problems of a surface crack in a structural 
component which may locally be represented by “plate” or “shell” have long 
attracted the attention of investigators who have been interested in fracture 
mechanics. In the past few years there has been some renewed interest in the 
Line-spring model, (which was originally developed by Rice [l]) for solving 
surface crack problems in plates and shells. This interest seems to be justified 
because of the relatively high accuracy of the results obtained from this model 
(see, for example, [2], [3] and [4] ). The purpose of this study is to extend the 
method to orthotropic elastic plates. The stress intensity factors for two 
orthotropic plates and one isotropic plate are extensively studied for through 
and part-through cracks. Here the isotropic case has been considered for the 
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purpose of comparison. In materials undergoing fracture it is often observed 
that small pre-existing cracks or cracks generated during service join together 
and grow into major cracks to cause final failure. Thus, theoretical analysis of 
the interaction of multiple cracks in the structure is quite important. 
Consequently, the emphasis in this study will be on the interaction of two 
collinear, semi-elliptic surface cracks or through cracks in an infinite plate 
subjected to uniform bending or tension. This problem is formulated by using 
Reissner’s plate bending theory and considering the material orthotropy with the 
isotropic medium as a special case. The problem is reduced to a system of 
singular integral equations. The effect of material orthotropy on the stress 
intensity factor and the interaction between two arbitrarily arranged collinear 
cracks have been investigated. Rather extensive numerical results have been 
provided which may be useful in application. 


Chapter 2 

The general formulation of the problem 


2.1 The governing equations for an orthotropic plate under 
bending and tension 

In this section a brief derivation of the fundamental equations for an 
orthotropic plate under bending and in-plane stretching is given. The problem of 
bending and the problem of in-plane loading are uncoupled. The material of the 
plate is assumed to be orthotropic with principal axes or orthotropy parallel to 
the co-ordinate axes. Thus the strain-stress relations may be expressed as 


€ = „ (J 4 S nn (T 4- 54^(7 

yy 21 zz 22 yy 23 zz ’ 


€ ~ 4- 4- S~ 0 <r , 

zz 31 ii 32 yy 33 zz 1 


e = S AA a , 

yz 44 yz ’ 


(i) 


C = S-r<T , 
zz 55 zz ’ 


zy 66 zy ’ 


where the quantities characterize elastic properties of the plate material. 
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2.1.1 The plate under bending 

The bending problem is based on the Reissner’s plate theory [5], [6]. In 
the isotropic case the equations are given by [7]. Consider a thin elastic plate 
bounded by the planes z=±h/2. The resultant forces are shown in Fig. 1, 
where M~ and V. (i,/=z,y) are respectively the moment and transverse shear 
resultants. 

Referring to [8], the equilibrium equations may be expressed as 


dM dM 

x x y 

+ -V- = 0, 


dx 


dy 


dM dM 
y *v .. 

dy + dx ~ ’ 


dV dV 
* v 

dx + dy ^ 


g = 0, 


( 2 ) 


and, in terms of the moment and shear resultants, the stress components 
are given by 


M z 
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a = — 
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v 
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M z 
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3V y 4z 2 
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3 q 2z 8z 3 2 
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(3) 


Also, the following relations between some average values w, 4> and <f> 

x y 

and the displacements u Q , v Q and are introduced in accordance with the 
balance of the work done by the resultant forces on the average values and by 
the corresponding stresses on the actual displacements 
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3 t +h ! 2 4 z 2 

“ = TkLn 


6 f +h > 2 z 


( 4 ) 


6 f +h / 2 z 

^y^^J-h/2 v °h/2 dz ' 


From (1) and the strain-displacement relations it follows that: 


d«i 0 

f X* = aT = S u °zz + S n a yy + S 13°zz > 

dv o 

i w = ~fy = S 21 ff « + S ™°W + S 23°zz > 

dw 0 

l zz = ~d 7 = S 31°xx + S 32°yy + 5 33 0 zz ’ 


du Q 

€ « = i>7 + d7 _S55<r «’ 


( 5 ) 


a« 0 <^ 0 _ 

( *V = ~dy ^Tx~ 5 66 °zy ’ 


V~ dy + dz ~ 44<T «*' 


If we now solve (5) for a^, ( t, j = i, y, z ), substitute into (3), integrate in 2 
and use (4) we obtain: 
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h , ™ z 6 <J 
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dx 


5h' 


h 3 , d * x d *y 

M = — (5,J _1 — + — ), 

xy 12 V 66' ' dy dx ’ 

dw 6V x 

** = -di +s ™Jh' 


( 6 ) 


du 6V y 

~ - dy + 5 44 5 h ' 


Technically, equations (2) and (6) would give the complete formulation of 
the problem, the eight equations accounting for the eight variables, M , M , 

x y 

M , V , V , <t> , <t> and u\ Eliminating first, 6 and 6 then M , M and M 

xy ’ x ’ y’ ^x’ ^ y 6 ’ r x r y 2’ y xy 

and assuming q— 0 , we obtain the system of equations governing the small 

deflection of elastic orthotropic plates as follows: 

V +V =0, 
v,y ’ 


V 4 -A~V ”h — A a u) — 0 , 

x 1 x,xx 2 x,yy 3 ,xxx 4 ,yyj 7 


( 7 ) 


V+A,V + A,V - A-U> — Arjjj =0, 

y 5 y,xx 6 y,yy 7 ,xxy 8 ,yyy 


where the coefficients A. are constants defined in Appendix 1. The 
remaining unknowns can be obtained from 
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55 


** = - w ,* + l T v *’ 


65 44 

+ 5/1 V 


M z = ~TT k ( 5 22 W ,xz " 5 12 W >W > + ( 5 22 5 55 + S 12 5 44 ) V *,z ’ 


M y = “ n~K ^ S J 1 W .SW " S 12 W ,« ) + 10 k ^ S 11 S 44 + S 12 5 55 ) V y,y > 


M zy $S 66 W ' X V + 10S 66 ^ S M V *’V + 5 44^y,i ) ’ 


(8) 


where k = S n S 22 ~ ^12^12 ' 

Note that since the basic system of equation (7) is sixth order, three 
conditions must be prescribed on the boundary of the plate. 

Following |8], to simplify this bending problem a stress function F(z,y) is 
introduced as follows: 


W = F ,y +A l F ,zxy +A 2 F ,yyy' 

^ 1 — ^3^,xxxy ^4^,xyyy ’ 


V y ~ *3 F ,xxxx ^4 F ,xxyy' 


(JO) 


then, equations (7) are identically satisfied provided the stress function 
F(z,y) satisfies the following equation: 
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( 11 ) 


d 4 F d 4 F d 4 F d 6 F 
B . — - + B„ — ~ — - + B , — + B. - 

1 dx 4 dx 2 dy 2 dy 4 dx* 


d*F 


d*F 


4 Br — “ — “ 4 JEh, ~ j 4 B 7 ~~~ = 0 j 

dx 4 dy 2 6 dx 2 dy 4 *- # 


d 6 F 

d y 6 


where Bj are known constants defined in Appendix 1. 

Equation (11) is called the fundamental equation. It is essentially a special 
case of that given in [8]. 


2.1.2 In-plane stretching problem 

Begining with the usual plane stress assumption 

a —a —a = 0 , 
zz zx zy ’ 

we define the Airy stress function $>(x,j/) as follows: 


( 12 ) 


d 2 $ 


d 2 * 


d l $ 


a XX ~ 2 

dy 




a = - 


x « dydx 

Along with the stress-strain relation (1), the compatibility condition for an 
orthotropic plate may be expressed as: 

d 4 $ 


d 4 $ „ d 4 * 

7+2C ~ “+C 2 7=0, 

dx 4 1 dx 2 dy 2 1 dy 4 


(13) 


where 


S 66 5 12 


S 1/2 


vv A* ^ * A 

C i = 2 5Z + sZ ’ C 2 = ( T ~ ) 


2S 22 b 22 


22 


(14) 


The fourth order differential equation (13) is subject to two conditions on 
the plate boundary. 

As we know, the solutions satisfying equation (13) yield the following 
characteristic equation 

m 4 4 2C 1 m 2 4C 2 2 = 0. (13.1) 
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The roots of (13.1) are 
™ 1 = ? 1 + *< 7 2 

"*2 = 93 + * </ 4 ( 13 2 ) 

m 3 — -m 1 , — ~ m 2 

Where m 1 and m 2 are both real or complex conjugates. The material is 
defined as type 1 when both roots are real and as type 2 when they are 
complex conjugates. 

Equations (11) and (13) give the governing equations for an orthotropic 
plate subjected to bending and in-plane streching. 

2.2 Infinite plate containing a series of collinear cracks 

2.2.1 The orthotropic plate under bending 

Consider now an infinite plate containing a series of collinear through 
cracks which are along the y axis. This problem will be solved by using Fourier 
transforms. It will be assumed throughout this study that through a proper 
superposition, the original crack problem has been reduced to a perturbation 
problem in which self-equilibrating force and moment resultants acting on the 
crack surfaces are the only nonzero external loads. One may note that since the 
corresponding through crack problem is uncoupled through the governing 
equations as well as the boundary conditions, the in-plane loading and bending 
problems can be solved independently. 

For the function F(x,y) defined by (11) we introduce the following Fourier 
transformations: 
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(15a) , 


/ +OO 

F(x,y) t iay dy , 

-oo 


1 f + °° 

F ( I,y ) = 2^7-00 ^ X,a ^ e ,aVdy 


(156) 


Substituting (15b) into (11) one obtains a sixth-order ordinary linear 
differential equation in f. Looking for a solution in the form: 

f(x,a) = R(a)e mz , (16) 


the charcteristic equation of the problem is found to be 
S 4 m 6 + ( flj - a 2 B 5 ) m 4 - ( a 2 B 2 - a 4 B g ) m 2 + ( a 4 B 3 - a*B 7 ) = 0 . 


(17) 


It should be emphasized here that the roots of (17) are in general complex 
and are functions of a. 

After solving (17) let the roots be ordered such that 
R t (mjl <0, (j=l, 2,3), m j+3 =- mj . (18) 


The solution /(*,<*), satisfying the regularity condition at x = ± oo, may 
then be expressed as : 

3 


f(x,«) = '£ t R j (a)e m f, for (*>0) 


;=i 


= £ for (l<0) - 
J=4 


(19) 


Because of the assumed symmetry with respect to yz plane in loading and 
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geometry, the transverse shear and moment resultants must satisfy the following 
symmetry conditions: 


v x (*,y) = ^,(*>1/) = V^(-x.y) , 

M x (x,y) = M z (-x,y), M y {x,y) = M y (~x,y ) , 


M xy (x,y) = -M xy {-x,y) . 


(20) 


It is therefore sufficient to consider one half of the plate only. Thus for 
example for x > 0 we find 

1 f°° 3 

F{x ' y) = uJ co £ ( * > 0 ) . (21) 

°°i=i 


Substituting from (21) into (10) and using (8) the relevant components of 
moment and transverse shear resultants and rotation may be expressed as 
follows: 


1 r°° 3 

M = — / Vr ( a)R{a)e m fc~ tyQ da , 
2nJ-oo Z ^ xyK 

3 = 1 


(22) 


9% ft Q ^ o A 

r “ I -i- / ■' A — C A ^ ( m ^ _L C A tr> ^ 1 /w * 

*V V • '6S 66 v v 65 66 1 105 66 “ 55 ' 105 66 "^" 4 V J “ * 

ft 3 A 2 h 2 

+ j — — — j4, -+ — — — 5,. ^4 . ] m -a 1 - — — — S. .4, ( m ) 5 , 

1 6 S cc 2 105^ 55 4) j 10S 44 3' ’ 
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( 23 ) 


i3 • r oo 3 

M = / Vr (a)J?.(a)e m /c-*‘' Q rfa, 

1 12/c 2nJ-oo xx ] 

;=i 


r„(») = I s 2! (rn/ + * 22 'l, (<■>/ - ^ (S 22 S S5 + S 12 S 4 < )^ 3 (m/|a 

6 

- [5 22 i4 2 (m j .) 2 -5 12 - 5 n j4j (m^) 2 - — (S 2 2 S 55 
+ 5 12 S 44 ) A 4 ( m y ) 2 ]a 3 - S 12 A 2 a 5 , 


r -\-OC ^ 

J=1 


( 24 ) 


*»~rJ! £ r » • 




( 25 ) 


6 ^55 


6 5 55 


r ,» = [»« i +Mi-g-jj-^3)»» i 1“ IV i"V* 4 > 


where Rj, R 2 and R 3 are unknown functions of a and are given in 
Appendix 2. 

In addition to the regularity condition at infinity, the bending problem 
must be solved under the following boundary conditions: 
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M zy (x,0) = 0, 

7 y (x,0)^0, -oc < x < oo, (26) 

A/ a (o f y) = p 1 (y), yei, 

<t> x { o,y) = o, y€L # , (27) 


where (L+Z/) = (-oo,+oo), L refers to a system of collinear cracks and is 
finite and Pj(y) is a known function . By using the homogeneous conditions 
(26), two of the three unknown functions R. can be eliminated. The third one is 
then obtained from the mixed boundary conditions (27). 

The problem may be reduced to an integral equation for a new r unknown 
function Gj defined by 


^<* 2 (0,y) = G^y), -oo <y< +oo. 


(28) 


Thus, from(27) by noting that Gj(y)=0 on y€ L' and by using (23) and 
(25) the first equation of (27) may be expressed in the following form: 


lim f G 1 (t)dtf + H(a,x)e tQ ^da= Pl (y), ( y eL). 

e — > 0 J L J -oo 


(29) 


The possible singularity of the kernels in (29) at y~t would be due to the 
behavior of H(a,x) as a — ► ±oo. For the purpose of examining the singular 
behavior of the kernels in (29) and for extracting the singular parts, the 
asymptotic behavior of rm, as |a|— » oo is needed. Thus, from (17) it can be 
shown that for large values of |a[ we have 
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m 1 fts a 1 |a|, 


m 2 « (a 2 + tfc 2 )|a|, m 5 = -m 2 , 

m 3 w ( a 2 “ 2 ) | a | , m 6 = -m 3 , 

(30) 

where a. and b. are material constants. 

Using now the relations (30) and separating the asymptotic value of H for 
large |a|, the kernel in (29) may be expressed as 


/ + oo r + oo 

H(a,x)e^-^ a da= / H °°( 0 , 1 ) 

-OO ^ “OO 


■*>da 


+ 



H°°(a,x) } e ia (‘-^da , 


(31) 


where H°° is the asymptotic value of H(a,x) for |q| — ► oc. The first term 
on the right-hand side of (31) gives a Canchy-type kernel (t-yp 1 and the 
second integral is uniformly converged for all t and y in which the limit x-0 
can, therefore, be put under the integral sign. After some rather complicated 
and lengthy manipulations (mostly due to the large number of elastic constants) 
(29) may be shown to be a singular integral equation with a simple Chachy- 
type singularity of the following form : 


fij |— + k{y,t)} G^tjdt = pj(y), ye L, 
J L r y 


( 32 ) 


where /ij is a material constant and 
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(33) 



H{a) - //°°(a)] S'l'-rida . 


If c % < y < d . , (t=l, . . . ,n) defines the cracks along the y axis, from 

(27) and (28) it follows that (34) must be solved under the following single- 
valuedness conditions: 

f 

fG (t)dt = 0. (t=l, . . . ,n) (34) 

Jc i 

For example, if the plate contains a single crack along (-a, a) subjected to 
uniform bending moment 

P,(») - -M°°, (35) 


it is seen that L=(—a,a) 
expressed as [10] 


and the solution of the integral equation may be 


*M0 

c,(0 = ^=- 

Vi-t 2 


(36) 


where the unknown function gj(t) is bounded in |t| < a. 
Defining now 
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T~t/a 


* — y/ a i 


h 2 ff b 


9l( T ) = 


h l( r ) 

VT?' 


<r L = 


6M°° 

U~’ 


equation (32) may be expressed as 


/: 


1 ?j( f ) 


Vl — T 2 


T — 8 


+ *( *,r ) ) dr= -1, |«j < 1, 


subjected to 


/; 


1 9i( T ) 


dr — 0. 




( 37 ) 


(38) 


(39) 


In the symmetric bending problem under consideration, the stress intensity 
factor at a crack tip y=<f. is defined as follows: 

*j(z)= lim \2(y- <r x (Q,y,z). (40) 

Refering to [10] for the procedure, in this single crack case the stress 
intensity factor ratio defined by 

*i6=*i(l/2 )/ k M = -j^ k bb a b y f< 1 ' ( 41 ) 


is found to be 
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( 42 ) 


We must emphasize again that the major difficulty in this problem is that 
the functions m^a) are not known explicitly in terms of (a). Thus we must 
make use of a numerical technique to obtain H( a) and then, k ( y , t ). 

2.2.2 The Orthopropic Plate under In-Plane Loading 
Similarly for the in-plane loading problem we let 

$(*>!/) = “/ tay da, (43) 

171 J-oo 


Substituting now from (43) into (13) and considering the symmetry 
conditions, we obtain: 

for materials of type 1 

<t>{x,a) = fijfa) exp(g 3 |a| * ) + « 2 ( a) exp(g 4 |a| x) y (x > 0) (44.1) 

where 

Wl l = »3 = -( C l + v/C l 2 - <7 2 2 ) 1/2 > 

™ 2 = <i 4 = -(c 1 -y/c 1 2 -c 2 2 ) 1 / 2 . (44.2) 


are real roots of the characteristic equation (13.1) and Cj and C 2 are 
given by (14). 

for materials of type 2 

<f>(x,a) = [ nj(a)cos ( q 2 ax ) + n 2 (a)sin ( q 2 ax ) ] expfaj |a| x ) , 

(x > 0) (45.1) 


where 
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m 


! = ?! + »' 9 2 < °) 

m 2 = — 1 


(45.2) 


are complex conjugate roots of the characteristic equation (13.1). 

The membrane problem must be solved under the following boundary 
condition: 


N xy (0,y) = 0, 


lim JV^O.y) = p 2 (y ) , 

i^0 + 


u(0 + ,y) = 0, 


ye L, 


y€L'. 


(46) 


(47) 


One of the two unknown functions and H 2 may be eliminated by using 
(46) and the remaining one may then be obtained from the mixed boundary 
conditions (47). Defining a new unknown function G 2 by 
d 


dy 


u(0 + ,i/) = G 2 (y), 


(48) 


using the following expressions obtained from (13), (43) and (44) or (45), 
for materials of type 1 
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( 49 . 1 ) 


r 

*V 


/ a[n i (a)g 3 |a| exp(g 3 |a|i) + n 2 (a)g 4 |a| exp(g 4 |«|i)l< tay da , 

Z7t J — oo 


a 

x 



-a 2 ) (nj(a)exp(g 3 |tt| 


*) 


+ fl 2 (a)exp(g 4 |a| x)}e lay da , 


for materials of type 2 




t 

= — / a { ^j( a ) [ — g 2 a sin (g 2 a * ) + 9j| ft | cos (g 2 a z ) ] 

Z 7T */ — oo 

+ fi 2 (a) | ? 2 a cos (g 2 a x ) + g^aj sin (g 2 a x ) ] } exp^ |a| x ) e~ ,ay da , 


( 49 . 2 ) 


1 f + °° 

0x= uJ-oo ^ _0,2 )I n i(°) COS ^2 al ) 

+ n 2 (a)sin (g 2 a x ) ] exp(g 4 |a| x ) e~ tay da , 


and noting that 


dr do do 

o u xy x y 

r = Sec “ — - ( 5. 9 — h Srffy ) J 

^ y2 <9 y ^2 d X 


( 50 ) 


one may obtain the following integral equation to determine G 2 


20 


f G 2<*> 

"4 T7'" = '’2 (,) ’ 


(y€ L) 


(51) 


where for materials of type 1 
h 1 

^ ' S 22^M + 

for materials of type 2 
h 1 

_ * 2 S 22 0 j ( 9j 2 + 9 2 2 ) 


(52.1) 


(52.2) 


Equation (51) must be solved under the following single-valuedness 
conditions, 



(• = 1 , 2 , 


(53) 


2.3 The line-spring model for surface crack 

The most general description of the line-spring model introduced in [l] has 
been given in [ll]. In a plate containing a part-through crack and subjected to 
membrane and bending loads, the net ligament (the uncracked portion) around 
the crack would generally have a constraining effect on the crack surface 
displacement. By approximately representing the net-ligament stress as a 
membrane load N and a bending moment M and the crack surface displacement 
due to N and M as an opening 6 and a rotation 0, the three-dimensional crack 
problem could be reduced to a two-dimensional coupled bending-membrane plate 
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problem. Furthermore, assuming that the relationship between (N,M) and (<S,0) 
may be found through the plane strain results obtained from the solution of an 
edge-notched strip, the pair of functions (£,0) and (M,N), then, are determined 
from the corresponding mixed boundary value problem for a plate having a 
through crack in which N and M are treated as unknown crack surface loads. 

Let the stress intensity factor for the plane strain problem in Fig. (2.b) be 
given as 

K(s) = y/h [ a + tn ] , (54) 

*{y) = L(y)/h, (55) 


where functious g b (s) and g f («) are called the shape functions for tension 
and bending, respectively and given as: 


g t ( 6 ) = y/xs^2 AT { * 2(t 1) , 
i=l 

n 

g b {s) = AB f 

t = l 


(56. o) 


(56.6) 


The coefficients AT. and AEk can be found by applying a suitable curve 
fitting to the results obtained from [12] which are valid for 
0 < Ljh < 0.8 . 


Referring to [13], the strain energy release rate Gj in an orthotropic 
medium may be obtained as: 
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1 1 1 
, * 11*22 , 2 , *22 2 2 * 12 + *66 2 
g, = *.’(-?-) Kr-) +-=[— I */<-3 


'11 


( 57 ) 


11 


where Kj is the stress intensity factor and the elastic constants b { . are 
givens by (see equation(l)) 




^22 ~ 


^ 33 S 22 ^23 


fc 66 _ 


S 22 

J 

S U S 22~ 

c 2 

b n 

S 22 


5 55 ’ 


S \3 S 22~ 

S 12 S 23 
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1 1 1 
b U b 22 ~ 2 b 22 2 2fe 12 + fc 66 ~2 


(58) 


Also, the energy available for fracture can be expressed as: 

d 1 dS 39 

C,« 5Z (C-V) = i (lV- + A,->. 


(59) 


1 

*3 


Combining now (57) and (59) and using (54) we obtain 

Id d fc 2 0 


( <r 2 S( 2 + 2<rm ?i ?, + m 2 g* ) = - | a— ( hS ) + m— ( — ) ] . 


(60) 


From (60) it can be shown that 


23 


a | a- 


( 61 ) 


e = 2h — ( a bb m + a H c ) , 

”3 


1 

6 = 2h ~ (a tb m+a u c), 
^3 


where a ^ are the compliance coefficients and are given by 


1 [ L 

i) = hJo 9 i 9 i dL ’ 


(62) 


From (28) and (48) it may be seen that 


0 — 2<£ 2 ( 0,3/ ) = 2 



<5 = 2u( 0,y ) = 2 / G 2 (t)d* . 




(63) 


Solving (61) for <r and m and substituting from (63), we obtain 


h 2 h 2 [» h [* 

A^(y) = — rr* = m 3 | — 7 w y G,(t)dt+- 1w ^ G 2 (<)^1, 

t I 

h fV f y 

^ (y) ~ h<r — fi 3 \ — 7 tb J' <?,(<)* + 7 tt c 2 («) rft l> 


(64) 
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1 bb - A > 7 M ~ A ’ 

“tfc _^U 

"Itb - ~ A ’ - A ’ 

A = a W a tt~ Q M 2 - 


( 65 ) 


One should keep in mind that consequently 7 ^ are fuctions of L(y)/y 
and that L(y) describes the surface crack profile. 

As an example let us now consider a plate containing a single surface 
crack, as shown in Fig. 2. By using this model, the crack may be assumed to 
be a through crack of length 2a and the constraint caused by the net ligament 
may be accounted for by applying the membrane and bending resultants N(y) 
and M(y) on the crack surfaces. One may note that these net ligament stresses 
tend to prevent the crack face from opening and rotating. The mixed boundary 
condition corresponding to (27) and (47) may be expressed as : 

M xx = G i(y) + ’ -a < y < a 

( 66 ) 

<^(0,y) = 0, |y| > a 


N XX (M = G 2 M + 1 -a < y < a (67) 

u(0,y) = 0. |y| > fl 


The integral equstions, corresponding to (32) and (51) then become : 
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Ivl < « 


(68. a) 


Mi J (^ + *(»><)! G i( < )* = p 1 (v) + -W(»). 

M 2 J \ — \G 2 (t)dt = p 2 (y) + N(y), |»| < « (68.6) 


where 


M(y) = M 3 I M lbb J- a Gl ^ dt + 6 7 * L a G 2^ dt ^ ( 68c ) 

*(») = M 3 I J 7« /_" G, (t) dt + 7 „ /_" G 2 (t) dt } . (68. J) 

As seen from (68) in the surface crack problem the integral equation for 
bending and tension are coupled. 

2.4 Two special cases of surface cracks 

As mentioned earlier, integral equations for the general case given in (32) 
and (51) are valid for any number of collinear cracks. But in this study, for the 
purpose of simplification, we only consider two special cases: a single symmetric 
crack and two arbitrarily arranged collinear cracks. 

In the first case, referring to Fig. 2, assuming the plate subjected to the 
uniform bending M°° and tension N°° and introducing the following change in 
variables 
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r = - , —a < (t,y) < a , 

a 

»= - , -1 < (r,«) < 1 , 

a 

G <(0 = ?,(’■) . (* = 1 , 2 ) , ( 69 ) 
* ( v,< ) = M «,0 , 


from(68), we obtain: 


A 2 /■* V 1 1 

J_ i 9 1 (r)dt- J _J — + afe,(«,r)] ?1 (r)dr 


h 

+ ^36 7 M a 



r= A/ 00 , 


(70.a) 


h f f 1 1 /•* 

J_ 1 9 2 {r)dT=N 00 . 


(70.6) 


It is seen that (70) has a simple Cauchy kernel and hence the solution of 
the integral equation is of the following form : 

M r ) 

?,(*■) =-=, (* = 1,2), -l<r<l. (71) 

Vl - r 2 


Equations (70) may be solved numerically by using the quadrature 
formulas given , for example , in [9], under the additional conditions: 



(* = 1 , 2 ). 


(72) 
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After solving (70) for g x and g v the stress and moment resultants M(y), 
N(y) may be obtained from (64), and the stress intensity factor K(y) along the 
crack front is then determined from (54). 

The configuration and parameters of the second special case are shown in 
Fig. 3. In this case the unknown functions G. should be defined separately for 
cracks (A,B) and (C,D) as: 

G n : derivative of the crack surface rotation for the crack (A,B); 

G 21 : derivative of the crack opening displacement for the crack (A,B); 

G 12 : derivative of the crack suface rotation for the crack (C,D); 

G 22 : derivative of the crack opening displacement for the crack (C,D); 
also let N(x) and M(x) for these two cracks be defined separately as follows : 
Nj(x) : N(x) for the crack (A,B); 

N 2 (x) : N(x) for the crack (C,D); 

Mj(x) : M(x) for the crack (A,B); 

M 2 (x) : M(x) for the crack (C,D). 

For crack (A,B) the system equation can then be written as: 

f B 1 

/ lrT" + MM)] G n( f )* 

J A 1 

- "i X lrr + 1 1 »!•' ) I c nW = 

*,(»,) - <■, J‘\ ^ ] C ,1<0* - n fc 71^ a nW = (") 

A < y l < B 
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similarly for crack (C,D) we have 


A/ 2 (y 2 )-/*!^ I^j^+ * ( 1 G n^ dt 

f D 1 

~Wc 


N 2 (y 2 )-fi 2 





- -G 22 (t)dt = N~, 
1 y 2 


C <y 2 < D 


(74) 


where 


h? h f 

•Mj (y) = I ~ i bb J' G u (0 + g i bt J G 2 i (0 dt I ’ 

N i (y) = m 3 1 1 7 (6 / Vt G lt (0 ^ V * G 2i (0 ] (75) 


Again, we have to convert the limits of integrals to (-1,1) for mumerical 
solution by introducing the following parameters: 
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B-A 


*1 = 

2 f * + 

2 =arj + (a + d/2), 


B-A 

B+A 

*1 = 

2 + 

2 — o*j + ( o + d/2 ) , 


D-C 

0+C 

t 2 = 

~2~ T * + 

2 -cr 2 + (c + d/2) , 


D-C 

ZHC 

* h = 

-6 ? + 
2 2 

2 - « 2 + ( c + rf/2 ) , 


= 9i/ r fc)> 

k 4 T i> e i) = 

TT ta -(y) = A w («) . (i,j, *=1,2) (76) 

Substituting from (75) into (73) and (74) we find 
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1 


h 2 /' 1 f l 1 

a#i 3 m /_ 1 l7^ + a *n ( r i-«i) 1 *n ( 7 i ) d7 i 

~ 1 ^ 72,41 ^ 9n ^ T? ^ ^ 7 2 + 0/i3 6 ? 2l( 7 l ^ ^ 7j = ^°° ’ ( 77 ) 


-1 < «j < 1 


h fl f 1 1 

- I 1 ? j ^2( r 2’®l) ^22^ 7 2^ T 2 + 0,l 3 X ttj_ } ? 2l( 7 l^ 7 l = ^°° > 


— !<«,<! 


(78) 


/l 2 f*2 f 1 1 

C/i 3 36 * bb J - 1 9 12^2^ 7 2 “ **1 y_jl r _ j c ^22^ 7 2’* 2^1 ^ 1 2 1 r 2 ) ^ 


f h f 2 

?22( 7 2 ) <f7 2 =:A/TO > 


(79) 


-1 < « 2 < 1 


(• I < 


^3 6 V i: ffl 2 ( 7 2) rf7 2 - ^2 1 ?22(7 2 )<f 


#* 2 f_! ^4^ 7 1’*2^ ff 2l( 7 l^ 7 l + C,l 3 \ it J } ^22^ 7 2^ 7 2 _ ^°° ’ 


(80) 


-1 < * 2 < 1 


where the kernels Fj, F 2 , F 3 and F 4 are bounded in the close interval 
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-l < ( V, ) ^ - 1 ’ 


(* = 1 , 2 ). 


This problem has four integral equations (77) -(80) with four unknown 
functions g yy g lT g 2l and g 22 which may be solved by using Gass-Chebyshev 
closed-type quadrature formula under the following single-valuedness conditions : 



(*,y=l,2). 


(81) 


After solving for the unknown functions, the same procedure for single 
crack may be used for each crack separately in order to find stress intensity 
factors. It may be observed that if we let AT = 0 , AT = 0 or 7 -= 0 , the 

integral equations (77)-(80) reduce to the system for two collinear through 

cracks in which the pair of equations corresponding to bending and tension 
would be uncoupled. One may also observe that as the distance between the 
two crack tends to infinity, the Fredholm kernels Fp • • • ,F 4 vanish and the 

integral equations corresponding to cracks (A,B) and (C,D) would again be 

uncoupled. For through cracks in the latter case one would obtain four 
uncoupled equations. 

2.5 Some remarks on the formulation of the orthotropic plate 
problem 

In the previous section the general problem for an orthotropic plate is 
formulated. The formulation is given for an orthotropic material which is 
defined by Hook’s law as expressed in (l). Since this is the first attempt to 
solve this kind problem, it is important to make the following remarks: 
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1. Unlike the classical orthotropic plate theory, here six rather than four 

material constants have been used to formulate the plate bending problem: four 

to describe the in-plane deformations names E , E , v , and G , (in terms of 

x y xy xy 

engineering material constants) and two to describe the out-of-plane 
deformations, namely G xe and G y£ . In the classical bending theory G xi and G yi 
are assumed to be infinite. 

2. Making use of the line-spring model in addition to the plate problem, 

an approximate solution is given for the three dimensional part-through crack 
problem. To do this the remaining three material coustants E e , j/ rx , 

(referring to Fig. 2) have also been introduced to the formulation. Thus, all 
nine material constants of orthotropy must be specified in order to solve the 

surface crack problem. 

3. As we mentioned earlier, if the plate material is isotropic, instead of 

one sixth order fundamental equation given by (11), the governing equation may 
be expressed in terms of one fourth order and one second order differential 

equation (see, for example [2]). Thus, it is analytically impossible to arrive at 
the same isotropic plate formulation as given in [2] by taking the limiting case 
of the orthotropic plate formulation given in this study. However, good 
agreement is obtained with the result given in [2] by using elastic constants in 
the numerical solution of the orthotropic plate problem which are nearly the 
same as the isotropic constants. It may also be noted that if the isotropic 

material constants are substituted into the characteristic equation (17), as 
expected the roots turn out to be the same as those found in the isotropic case. 
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Chapter 3 

Results and discussion 

The main interest in this study is in evaluating the stress intensity factor 
in an orthotropic plate. The similar study for the isotropic plate has been 

considered before in [2] and [10]. However, some basic and additional results 
for this case have been studied for the purpose of comparison and extension. 

The elastic constants of the materials used in the numerical examples are 
given in Table 1. Material 1 is basically isotropic, and 2 and 3 are laminated 
composite material. Except for a 90 degree rotation of the axes of orthotropy, 
Material 2 and 3 are identical. Most of the results obtained are for loading 
conditions 

JV(y) = 7V°°, A/(y) = 0, (tension) 

M(y) — M°° y N(y)-0. ( bending ) 

The first case that is studied is the plate containing a through crack and 
two collinear through cracks. Because of the nature of the plate theory used in 
analysis, the stress intensity factor is a linear function of the thickness 
coordinate 2 [see (41)]. 

Table 8 and Fig. 4 show the effect of the thickness ratio on the stress 
intensity factor for a single through crack. The result for the isotropic case 
(i/ = 0.3) are the same as that in [10]. It is clearly seen that the orthotropic 
materials have relatively higher stress intensity factors . This is due to a higher 

in these orthotropic cases. 

Some numerical results obtained for two arbitrarily arranged collinear 
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through cracks in a plate under uniform bending are shown in Tables 9-14 and 
Figures 5-6. Tables 9-11 show the interaction between two cracks with the 

length of one crack varying for Material 1, 2 and 3 and Tables 12-14 show 
similar results with the plate thickness varying. Tables 15-16 show the results 
for a plate under uniform tension. As expected for the distance between the two 
cracks d — ► oo the single crack results are recovered, and the stress intensity 

factor at the inner tip is always greater than at the outer tip. Figures 5 and 6 

clearly show the effect of bending due to the nonsymmetric orientation of the 
two cracks in the y direction. It is again noted that for the in-plane stretching 
problem the material orthotropy has no effect on the stress intensity factor for 
an infinite orthotropic plate containing a single crack or series of collinear 
cracks. 

After through cracks, the part-through crack problem is considered. As 
noted before, for the application of the line-spring model, the contour of the 
part-through crack can be any reasonable curve. Elliptic cracks are studied here 

since it is believed that the ellipse is the closest contour for the actual shape of 

the crack which may be encountered in practical applications. Thus, crack 
length for any cross section , referring to Fig. 2(a), is defined by 

L (y) = L 0 Vl - (y/a ) 2 , —a<y<a, 

L q being the total crack depth at the midsection (y = 0). 

Tables 17-19 show the stress intensity factor at y=0 in an infinite plate 
containing a single crack under uniform tension and pure bending for Materials 
1, 2 and 3 respectively. In these cases, the stress intensity factors are 
normalized with respect to K°° which is the corresponding value for an edge 
crack strip under plane strain condition with the same LJh ratio [Fig. 2(b)]. 
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Note that the limiting values of these stress intensity factors are 


0 

1 

S*5 

for 

K — K °° , 

for 



0 , 



oo . 


It should be noted in this case that because of crack closure on the 
compression side, taken separately the bending results are meaningless. 
However, they may be used by superposition with tension results which are 
sufficiently large so that the stress intensity factor on both sides of the crack 
are positive. The coefficients of the shape functions g b and g % for these three 
materials are listed in Tables 3, 5 and 7 which is based on the stress intensity 
factors shown in Tables 2, 4 and 6. 

Figures 7 and 8 clearly show the comparison of different material on the 
stress intensity factor for Ljh varying from 0.2 - 0.8. As expected under 
bending the crack faces would be closed on the compression side; therefore, for 
deep cracks the stress intensity factor in an isotropic medium becomes negative. 
From Tables 18 and 19 and Figures 8 and 9 we have found that the stress 
intensity factor obtains the negative value as the LJh larger than 0.5. 

Some results regarding the distribution of the normalized stress intensity 
factor at the crack front for a single semi-elliptic surface crack are given in 
Table 20. 

The results for two collinear semi-elliptic surface cracks are shown in 
Tables 21-26. Here the crack profile is defined by 
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L(y) = L 0 V 1-y 2 , 


V - (y B + V A ) I 2 y-(B+A)/2 

v ~ (t/B-yJ/ 2 " (B-A)/2 


y-(yp + y c )/ 2 y -{D+c )/2 

y= (y/>-y c )/ 2 “ ’ c< y <D 

Results for the interaction of these two collinear, semi-elliptic surface 
cracks with the distance between them varying from d/a= 0.1 - oo and the 
length of one crack varying from e/a =0.1 - 1. under tension or bending have 
been tabulated in Tables 21-23 for material 1, 2 and 3 respectively. The effect 
of Ljh on the normalized stress intensity factor for the case of these two 
collinear cracks has also been considered. Tables 24-26 show this effect for the 
case of two identical collinear cracks. It may be observed again from all of 
these tables that the material orthotropy does have significant effect on the 
fracture parameter for the surface crack. 

Finally, an additional example has been given for a crystalline material 
TOPAZ [ Sio 4 AL 2 (F,OH) 2 ] which has an orthohomibiz system behavior. The 
material constants of TOPAZ, given by [14], are as follows: 

(unit 10’ 13 cm 2 /dyne ) 


S n =4.43 S 22 =3.53 S 33 =3.84 

Material 4 S 44 =9.25 S &5 =7.52 S 66= 7 ' 63 

S 12 =-1.38 S 13 =-0.86 S 23 — -0.66 
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Material 5 


S u =3.53 S 22 =4.43 S 33 =3.84 

S 44= 7 - 52 S 55= 9 - 25 S 66= 7 - 63 

S 12 =-1.38 S 13 =-0.66 S 23 =-0.86 

Note that the only difference between Material 4 and Material 5 is a 90 
degree rotation of orthotropy. Table 27 shows the effect of the thickness ratio 
ajh on the stress intensity factor in a single through cracked plate under 
uniform bending for Material 4 and Material 5. and Table 28 shows the effect 
of Ljh on the normalized stress intensity factor at the deepest penetration 

i 

point of a semi-elliptic surface crack under uniform tension and bending 
moment. It may be seen from these tables that there is no important difference 
between these two kinds of orthotropy orientations, as this material does not 
have strong orthotropic material properties. 

! 

i 

! 

| 

i 
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Tables and Figures 


Table 1: Elastic constants of the materials used in the examples 


Mat.l Mat.2 Mat. 3 


E x psi 

GPa 

2.2447-10 6 

15.477 

5.8610 6 

40.405 

22.2-10 6 

153.069 

E y psi 

2.2610 6 

22.2-10 6 

5.86-10 6 

GPa 

15.583 

153.069 

40.405 

E , Psi 

2.2610 6 

3.3-10 6 

3.310 6 

GPa 

15.583 

22.754 

22.754 

G xy P S > 

0.866-10 6 

4.25-10 6 

4.25-10 6 

GPa 

5.971 

29.304 

29.304 

G y* P si 

0.866 10 6 

0.592 10 6 

0.225-10 6 

GPa 

5.971 

4.082 

1.551 

G xz P si 

0.866 10 6 

0.225- 10 6 

0.592 10 6 

GPa 

5.971 

1.551 

4.082 

V 

0.3 

0.484 

1.834 

xy 

V 

0.3 

0.195 

0.261 

yz 


0.3 

0.261 

0.195 
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Table 2: The coefficients of AT ; and ABj for the shape 

functions g t and g b ( Material 1) 


AT. AB; 


1 

1.122 

1.120 

2 

6.520 

-1.887 

3 

-12.388 

18.014 

4 

89.055 

-87.385 

5 

-188.608 

241.912 

6 

207.387 

-319.940 

7 

-32.052 

168.011 


Table 3: 


Stress intensity factors in an edge crack under 
tension (N) or bending (M), (Material 1) 



6Af 

l= :r > 


L 

h 


<t^L 



nW~L 


0.001 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 


1.122 

1.189 

I. 367 
1.660 
2.111 
2.825 
4.033 
6.355 

II. 955 


1.120 

1.047 

1.055 

1.124 

1.261 

1.498 

1.915 

2.728 

4.691 
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Table 4: The coefficients of AT ; and AB. for the shape 

functions g t and g b ( Material 2) 


AT, AB, 


1 

1.047 

1.043 

2 

7.639 

-1.610 

3 

-27.969 

17.276 

4 

175.360 

-84.989 

5 

-439.451 

232.556 

6 

557.540 

-304.196 

7 

-222.80 

158.307 


Table 5: 


Stress intensity factors in an edge crack under 
tension (N) or bending (M), (Material 2) 


(<r= 


N 

I* 


6Af . 



L 

h 


K. 


K 


M 


o'Tl 


m^L 


0.001 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 


1.042 
1.129 
1.318 

I. 607 

2.042 
2.721 
3.860 
6.038 

II. 277 


1.041 

0.992 

1.013 

1.083 

1.211 

1.430 

1.814 

2.563 

4.372 
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Table 6: The coefficients of ATj and AB. for the shape 

functions g t and g b ( Material 3) 


AT. 


AB. 


1 

1.055 

1.051 

2 

7.461 

-1.664 

3 

-25.426 

17.517 

4 

160.659 

-85.711 

5 

-396.808 

234.182 

6 

498.577 

-306.252 

7 

-191.359 

159.403 


Table 7: 


Stress intensity factors in an edge crack under 
tension (N) or bending (M), (Material 3) 




L 

h 




m^L 


0.001 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 


1.050 

1.134 

I. 321 
1.611 
2.048 
2.730 
3.876 
6.068 

II. 350 


1.049 

0.996 

1.016 

1.086 

1.215 

1.436 

1.823 

2.577 

4.402 
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Table 8: The effect of the thickness ratio on the stress intensity 

factor in a plate under uniform bending, containing 
a single crack, a b = 6M Q /h 2 


Mat.l Mat. 2 Mat. 3 


0.05 

0.9885 

0.9992 

0.9987 

0.1 

0.9676 

0.9974 

0.9957 

0.25 

0.8992 

0.9886 

0.9823 

0.5 

0.8195 

0.9689 

0.9548 

1 

0.7477 

0.9276 

0.9050 

2 

0.7003 

0.8664 

0.8434 

4 

0.6701 

0.8086 

0.7942 

6 

0.6446 

0.7776 

0.7657 

10 

0.6481 

0.7166 

0.7044 



43 


Table 9: Stress intensity factors in a plate containing two collinear cracks 

subjected to uniform bending moment M Q ( Material 1) 

. B-A D-c a . 

( d=C-B, -=1 ) 



c 

a 

— > 0.1 
a 

0.25 

0.5 

1 

2 

OO 


l 

0.8799 

0.8551 

0.8313 

0.8045 

0.7798 

0.7477 

^1A 

0.5 

0.8071 

0.7938 

0.7821 

0.7698 

0.7593 

0.7477 



0.25 

0.7711 

0.7647 

0.7598 

0.7551 

0.7513 

0.7477 


0.1 

0.7532 

0.7512 

0.7500 

0.7490 

0.7482 

0.7477 


1 

1.294 

1.076 

0.9599 

0.8697 

0.8049 

0.7477 

k -[R 



0.5 

1.063 

0.9143 

0.8458 

0.7995 

0.7698 

0.7477 

a 


0.25 

0.9161 

0.8220 

0.7863 

0.7663 

0.7550 

0.7477 


0.1 

0.8088 

0.7678 

0.7563 

0.7514 

0.7489 

0.7477 


1 

1.294 

1.076 

0.9599 

0.8697 

0.8049 

0.7477 

k 1C 

0.5 

1.012 

0.8405 

0.7498 

0.6786 

0.6261 

0.5794 

ft 


0.25 

0.7990 

0.6595 

0.5867 

0.5297 

0.4872 

0.4496 


0.1 

0.5647 

0.4577 

0.4037 

0.3627 

0.3325 

0.3060 


1 

0.8799 

0.8551 

0.8313 

0.8045 

0.7798 

0.7477 

k \D 

0.5 

0.7395 

0.7071 

0.6771 

0.6434 

0.6132 

0.5794 

c^a 


0.25 

0.6275 

0.5867 

0.5507 

0.5135 

0.4816 

0.4496 


0.1 

0.4817 

0.4293 

0.3917 

0.3577 

0.3308 

0.3060 
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Table 10: 


: Stress intensity factors in a plate containing two collinear cracks 

subjected to uniform bending moment M Q ( Material 2) 


, B-A D-c , ^ ^ a 
( > c “ — ? ? d=C-B, ~ ~1 ) 



c 

a 

-^ 0.1 

a 

0.25 

0.5 

1 

2 

oo 


l 

1.051 

1.023 

1.001 

0.9782 

0.9590 

0.9276 


— 

0.5 

0.9814 

0.9658 

0.9546 

0.9451 

0.9373 

0.9276 



0.25 

0.9488 

0.9411 

0.9363 

0.9328 

0.9305 

0.9276 


0.1 

0.9329 

0.9305 

0.9293 

0.9285 

0.9280 

0.9276 



1 

1.591 

1.274 

1.122 

1.029 

0.9773 

0.9276 

*1 B 

— 

0.5 

1.320 

1.106 

1.015 

0.9665 

0.9441 

0.9276 



0.25 

1.137 

1.008 

0.9614 

0.9405 

0.9325 

0.9276 


0.1 

1.003 

0.9500 

0.9353 

0.9301 

0.9284 

0.9276 


*1C 

l 

0.5 

0.25 

0.1 

1.591 

1.221 

0.9261 

0.6327 

1.274 

0.9665 

0.7202 

0.4774 

1.122 

0.8424 

0.6188 

0.4023 

1.029 

0.7658 

0.5567 

0.3577 

0.9773 

0.7240 

0.5236 

0.3348 

0.9276 

0.6851 

0.4943 

0.3154 


1 

1.051 

1.023 

1.001 

0.9782 

0.9590 

0.9276 

* ID 



0.5 

0.8441 

0.8019 

0.7687 

0.7386 

0.7151 

0.6851 

q 


0.25 

0.6817 

0.6228 

0.5790 

0.5436 

0.5201 

0.4943 


0.1 

0.5143 

0.4375 

0.3880 

0.3535 

0.3335 

0.3154 
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Table 11: 


Stress intensity factors in a plate containing two collinear cracks 
subjected to uniform bending moment M Q ( Material 3) 


( 



d=C-B, 




c 

a 

--^ 0.1 

a 

0.25 

0.5 

1 

2 

OO 


i 

1.031 

1.004 

0.9821 

0.9582 

0.9375 

0.9050 

*14 

— 

0.5 

0.9595 

0.9442 

0.9332 

0.9238 

0.9152 

0.9050 

<r* a 


0.25 

0.9263 

0.9188 

0.9141 

0.9106 

0.9083 

0.9050 


0.1 

0.9194 

0.9080 

0.9068 

0.9060 

0.9059 

0.9050 



l 

1.547 

1.246 

1.101 

1.011 

0.9581 

0.9050 

*1 B 

— 

0.5 

1.287 

1.080 

0.9917 

0.9452 

0.9228 

0.9050 

<7 


0.25 

1.111 

0.9841 

0.9385 

0.9183 

0.9103 

0.9050 


0.1 

0.9103 

0.9271 

0.9127 

0.9076 

0.9054 

0.9050 


* 1C 

l 

0.5 

0.25 

0.1 

1.547 

1.206 

0.9205 

0.6308 

1.246 

0.9550 

0.7171 

0.4775 

1.101 

0.8301 

0.6150 

0.4019 

1.011 

0.7507 

0.5505 

0.3555 

0.9581 

0.7067 

0.5154 

0.3311 

0.9050 

0.6752 

0.4912 

0.3149 


1 

1.031 

1.004 

0.9821 

0.9582 

0.9375 

0.9050 

k lD 

— 

0.5 

0.8289 

0.7869 

0.7529 

0.7216 

0.6981 

0.6752 

aVfl 


0.25 

0.6777 

0.6187 

0.5736 

0.5364 

0.5116 

0.4912 


0.1 

0.5140 

0.4375 

0.3872 

0.3510 

0.3300 

0.3149 
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Table 12: 


The effect of the thickness ratio on the stress intensity factors 
under uniform bending moment M Q (Material 1) 


( 


B-A 


c=— — — , d=C-B, -=0.5 ) 

2 ’ a ’ 



a 

h 

-^0.1 

a 

0.25 

0.5 

1 

2 

OO 


0.5 

0.8715 

0.8596 

0.8502 

0.8408 

0.8320 

0.8195 

* IA 

1 

0.8071 

0.7938 

0.7821 

0.7698 

0.7593 

0.7477 

7?z 

2 

0.7625 

0.7463 

0.7322 

0.7191 

0.7094 

0.7003 


* IB 

a''/ a 

0.5 

1 

2 

1.148 

1.063 

1.039 

0.9757 

0.9143 

0.8898 

0.9031 

0.8458 

0.8071 

0.8631 

0.7995 

0.7493 

0.8409 

0.7698 

0.7184 

0.8195 

0.7477 

0.7003 


0.5 

1.082 

0.8825 

0.7872 

0.7256 

0.6841 

0.6358 

k lC 

i 

1.012 

0.8405 

0.7498 

0.6786 

0.6261 

0.5794 



2 

0.9896 

0.8147 

0.7061 

0.6204 

0.5675 

0.5287 


0.5 

0.7831 

0.7522 

0.7260 

0.6989 

0.6732 

0.6358 

*in 



1 

0.7395 

0.7071 

0.6771 

0.6434 

0.6132 

0.5794 

a 


2 

0.6933 

0.6549 

0.6192 

0.5837 

0.5561 

0.5287 
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Table IS: 


The effect of the thickness ratio on the stress intensity factors 
under uniform bending moment M Q ( Material 2) 



c=— — — , d=C-B, -=0.5 ) 

2 ’a ’ 



a 

h 

-^0.1 

a 

0.25 

0.5 

1 

2 

oo 


0.5 

1.027 

1.009 

0.9956 

0.9849 

0.9772 

0.9689 

k JA 

1 

0.9814 

0.9658 

0.9546 

0.9451 

0.9373 

0.9276 

a^a. 

2 

0.9226 

0.9087 

0.8958 

0.8890 

0.8865 

0.8664 


* IB 

cr'^a 

0.5 

1 

2 

1.390 

1.320 

1.227 

1.160 

1.106 

1.036 

1.061 

1.015 

0.9550 

1.008 

0.9665 

0.9100 

0.9836 

0.9441 

0.8771 

0.9689 

0.9276 

0.8664 


0.5 

1.279 

1.002 

0.8652 

0.7796 

0.7336 

0.6990 

k 1C 

1 

1.221 

0.9665 

0.8424 

0.7658 

0.7240 

0.6851 

o'J a 


2 

1.148 

0.9248 

0.8164 

0.7462 

0.7024 

0.6559 


0.5 

0.8692 

0.8214 

0.7834 

0.7495 

0.7250 

0.6990 

n 

— 

1 

0.8441 

0.8019 

0.7687 

0.7386 

0.7151 

0.6851 

a 

2 

0.8144 

0.7792 

0.7463 

0.7204 

0.6905 

0.6559 
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Table 14: 


The effect of the thickness ratio on the stress intensity factors 
under uniform bending moment M Q ( Material 3) 


( 


B-A 


a 


2 » 



d=C-B, -=0.5 ) 

a ' 



a 

h 

--0.1 

a 

0.25 

0.5 

1 

2 

oo 


0.5 

1.011 

0.9942 

0.9817 

0.9713 

0.9701 

0.9548 

k lA 

1 

0.9595 

0.9442 

0.9332 

0.9238 

0.9152 

0.9050 



2 

0.9006 

0.8867 

0.8720 

0.8671 

0.8530 

0.8433 


0.5 

1.371 

1.143 

1.046 

0.9937 

0.9638 

0.9548 

k IB 

1 

1.287 

1.080 

0.9917 

0.9452 

0.9228 

0.9050 

0 r^a 


2 

1.188 

1.009 

0.9323 

0.8880 

0.8638 

0.8433 


0.5 

1.274 

1.000 

0.8641 

0.7773 

0.7293 

0.6946 

k 1C 

1 

1.206 

0.9550 

0.8301 

0.7507 

0.7067 

0.6752 



2 

1.106 

0.8840 

0.7743 

0.7049 

0.6668 

0.6400 


0.5 

0.8671 

0.8193 

0.7809 

0.7460 

0.7200 

0.6946 

* ID 

1 

0.8289 

0.7869 

0.7529 

0.7216 

0.6981 

0.6752 

fc> 


2 

0.7688 

0.7341 

0.7056 

0.6791 

0.6592 

0.6400 
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Table 15 : 


Stress intensity factors in a plate containing two collinear cracks 
subjected to uniform tension N Q 


, B-A _D-C a . 

( a ~ o ? c j * d — C - B , - — 1 ) 



c 

a 

-- 0.1 

a 

0.25 

0.5 

1 

2 

OO 


i 

1.151 

1.112 

1.0811 

1.052 

1.028 

1.000 

k lA 

0.5 

1.066 

1.045 

1.030 

1.018 

1.009 

1.000 

<r y ^a 

0.25 

1.026 

1.016 

1.010 

1.005 

1.003 

1.000 


0.1 

1.007 

1.004 

1.002 

1.001 

1.000 

1.000 


* IB 

1 

0.5 

0.25 

0.1 

1.795 

1.449 

1.234 

1.083 

1.414 

1.206 

1.090 

1.025 

1.229 

1.100 

1.038 

1.009 

1.113 

1.043 

1.014 

1.003 

1.048 

1.016 

1.005 

1.001 

1.000 

1.000 

1.000 

1.000 


1 

1.795 

1.414 

1.229 

1.113 

1.048 

1.000 

k 1C 

0.5 

1.329 

1.035 

0.8889 

0.7962 

0.7446 

0.7071 

o y /~a 


0.25 

0.9915 

0.7595 

0.5955 

0.5683 

0.5280 

0.5000 


0.1 

0.6732 

0.5006 

0.4151 

0.3624 

0.3349 

0.3162 


1 

1.151 

1.112 

1.0811 

1.052 

1.028 

1.000 




0.5 

0.8956 

0.8429 

0.8007 

0.7626 

0.7347 

0.7071 



0.25 

0.7170 

0.6480 

0.6425 

0.5524 

0.5235 

0.5000 


0.1 

0.5422 

0.4557 

0.3985 

0.3573 

0.3336 

0.3162 
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Table 16 : The effect of the thickness ratio on the stress intensity factors 

under uniform tension N Q 

( o=— y-, C=££ d=C-B, ^=0.5 ) 



a 

h 

-- 0.1 

a 

0.25 

0.5 

1 

2 

OO 


0.5 

1.066 

1.045 

1.030 

1.018 

1.009 

1.000 

k lA 

1 

1.066 

1.045 

1.030 

1.018 

1.009 

1.000 


2 

1.066 

1.045 

1.030 

1.018 

1.009 

1.000 


k \B 

a^a 

0.5 

1 

2 

1.449 

1.449 

1.449 

1.206 

1.206 

1.206 

1.100 

1.100 

1.100 

1.043 

1.043 

1.043 

1.016 

1.016 

1.016 

1.000 

1.000 

1.000 


0.5 

1.329 

1.035 

0.8890 

0.7962 

0.7446 

0.7071 

Jfe 

K ic 

1 

1.329 

1.035 

0.8890 

0.7962 

0.7446 

0.7071 

<7^a 


2 

1.329 

1.035 

0.8890 

0.7962 

0.7446 

0.7071 


0.5 

0.896 

0.843 

0.801 

0.763 

0.735 

0.7071 

k ID 

1 

0.896 

0.843 

0.801 

0.763 

0.735 

0.7071 

ft 


2 

0.896 

0.843 

0.801 

0.763 

0.735 

0.7071 
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Table 17 : Stress intensity factors at the maximum penetration point of 

a semi*elliptic surface crack in an infinite plate under uniform 
tension or pure bending (Material 1) 


-- 0.2 0.4 0.6 0.8 

n 

a 

h 



0.1 

0.378 

K _u 

0.25 

0.585 

0.5 

0.715 

K°° 

1 

0.811 


2 

0.882 


0.060 

- 0.039 

- 0.0296 

0.184 

- 0.003 

- 0.0345 

0.309 

0.053 

- 0.0287 

0.446 

0.136 

- 0.0119 

0.586 

0.242 

0.0179 



0.1 

0.423 


0.25 

0.613 

K u 

— 

0.5 

0.733 

K°° 


1 

0.823 


2 

0.889 


0.176 

0.079 

0.021 

0.284 

0.124 

0.034 

0.392 

0.176 

0.050 

0.510 

0.247 

0.073 

0.631 

0.339 

0.104 
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Table 18: Stress intensity factors at the maximum penetration point of 

a semi-elliptic surface crack in an infinite plate under uniform 
tension or pure bending (Material 2) 


-0.2 0.4 0.6 0.8 

h 

a 

h 



0.25 

0.518 

0.130 

-0.025 

-0.036 


0.5 

0.675 

0.255 

0.019 

-0.036 


1 

0.798 

0.409 

0.097 

-0.026 

K n 

— 

2 

0.873 

0.559 

0.205 

-0.0005 

K°° 


4 

0.924 

0.693 

0.339 

0.045 


6 

0.944 

0.759 

0.425 

0.082 


10 

0.962 

0.829 

0.536 

0.142 



0.25 

0.552 


0.5 

0.697 


1 

0.810 

K u 

— 

2 

0.881 

K °o 


4 

0.928 


6 

0.947 


10 

0.965 


0.237 

0.102 

0.027 

0.345 

0.147 

0.039 

0.477 

0.215 

0.058 

0.608 

0.307 

0.086 

0.726 

0.422 

0.129 

0.784 

0.496 

0.163 

0.846 

0.591 

0.217 
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Table 10: 


Stress intensity factors at the maximum penetration point of 
a semi-elliptic surface crack in an infinite plate under uniform 
tension or pure bending (Material 3) 


—0.2 0.4 0.6 0.8 

h 

a 

h 



0.25 

0.629 

0.212 

0.0024 

-0.037 


0.5 

0.765 

0.361 

0.069 

-0.031 


1 

0.861 

0.521 

0.172 

-0.103 


— 

2 

0.913 

0.660 

0.299 

0.028 

K°° 


4 

0.948 

0.775 

0.446 

0.090 


6 

0.963 

0.829 

0.534 

0.139 


10 

0.975 

0.882 

0.640 

0.215 



0.25 

0.654 


0.5 

0.780 


1 

0.869 

2 

0.918 

K°° 

4 

0.952 


6 

0.965 


10 

0.977 


0.308 

0.131 

0.035 

0.435 

0.191 

0.051 

0.574 

0.279 

0.076 

0.696 

0.387 

0.114 

0.798 

0.513 

0.170 

0.846 

0.589 

0.214 

0.894 

0.681 

0.281 
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0.000 


intensity factors at the crack front 
iptic surface crack with a/h-1 


1 


Mat.l Mat. 2 Mat.3 


0.509 

0.498 

0.544 

0.481 

0.453 

0.538 

0.450 

0.416 

0.522 

0.446 

0.409 

0.521 


0.396 

0.391 

0.419 

0.453 

0.433 

0.498 

0.491 

0.461 

0.549 

0.510 

0.477 

0.574 


0.275 

0.254 

0.321 

0.211 

0.183 

0.257 

0.154 

0.121 

0.195 

0.136 

0.097 

0.172 

0.210 

0.200 

0.243 

0.224 

0.203 

0.258 

0.239 

0.212 

0.273 

0.247 

0.215 

0.279 
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Table 21: 


The interaction between two cracks on the normalized stress 
intensity factors at the deepst penetration point of a semi- 
elliptic surface crack under uniform bending or tension 

(Material 1) 


B-a d-c , ^ „ « 


VQab 


B+A 


^0 cd 


D+C 




c 

a 

— > 0.1 
a 

0.25 

0.5 

1 

2 

OO 


i 

0.313 

0.306 

0.299 

0.290 

0.283 

0.272 

*- ( B ) 

Qab 

0.5 

0.292 

0.287 

0.282 

0.278 

0.274 

0.272 

K °° 


0.25 

0.280 

0.275 

0.275 

0.273 

0.272 

0.272 


0.1 

0.273 

0.273 

0.272 

0.272 

0.272 

0.272 



i 

0.313 

0.301 

0.299 

0.290 

0.283 

0.272 

^oJ B) 

0.5 

0.197 

0.188 

0.179 

0.171 

0.164 

0.160 

K °° 

0.25 

0.101 

0.091 

0.083 

0.076 

0.072 

0.069 


0.1 

0.012 

0.0045 

- 0.0004 

- 0.0038 

- 0.0057 

- 0.0058 


1 

0.397 

0.392 

0.386 

0.379 

0.374 

0.366 

*0a6 (71 
• ,oo 

0.5 

0.382 

0.378 

0.375 

0.371 

0.368 

0.366 

A 

0.25 

0.373 

0.371 

0.369 

0.368 

0.366 

0.366 


0.1 

0.367 

0.367 

0.366 

0.366 

0.366 

0.366 

K 0cd (1) 

1 

0.397 

0.392 

0.386 

0.379 

0.374 

0.366 

0.5 

0.300 

0.293 

0.286 

0.279 

0.274 

0.269 

K °° 

0.25 

0.217 

0.209 

0.203 

0.198 

0.194 

0.190 


0.1 

0.136 

0.130 

0.126 

0.124 

0.124 

0.123 
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Table 22: The interaction between two cracks on the normalized stress 

intensity factor at the deepst penetration point of a semi- 
elliptic surface crack under uniform bending or tension 

(Material 2) 





B-A 

D-c , ~ 

a 





a 

— , c = 

, d = C 

- S , 7=1 






2 ’ 

2 ’ 

’ h 







B+A 

D+C L o 

= 0.4 





^0 ab 

= 2 ’ V 0cd 

2 ’ h 



e 

a 

--* 0.1 

a 

0.25 

0.5 

1 

2 

OO 

K ( B ) 

K 0ab 

i 

0.444 

0.437 

0.430 

0.423 

0.417 

0.409 

0.5 

0.430 

0.425 

0.420 

0.416 

0.412 

0.409 

K°° 

0.25 

0.419 

0.416 

0.414 

0.412 

0.410 

0.409 


0.1 

0.412 

0.411 

0.410 

0.410 

0.410 

0.409 

«cJ B) 

1 

0.444 

0.437 

0.430 

0.423 

0.417 

0.409 

0.5 

0.298 

0.287 

0.277 

0.268 

0.262 

0.255 

K°° 

0.25 

0.171 

0.158 

0.148 

0.140 

0.135 

0.130 


0.1 

0.054 

0.043 

0.036 

0.032 

0.029 

0.027 


1 

0.505 

0.500 

0.494 

0.488 

0.483 

0.477 


0.5 

0.494 

0.490 

0.486 

0.482 

0.479 

0.477 

K°° 

0.25 

0.485 

0.482 

0.480 

0.479 

0.478 

0.477 


0.1 

0.479 

0.478 

0.477 

0.477 

0.477 

0.477 

K (1) 
K Ocd 

1 

0.505 

0.500 

0.494 

0.488 

0.483 

0.477 

0.5 

0.380 

0.370 

0.362 

0.355 

0.350 

0.345 

K°° 

0.25 

0.270 

0.260 

0.252 

0.245 

0.241 

0.237 


0.1 

0.168 

0.160 

0.154 

0.150 

0.148 

0.146 
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Table 23 


: The interaction between two cracks on the normalized stress 

intensity factors at the deepst penetration point of a semi- 
elliptic surface crack under uniform bending or tension 

( Material 3) 


a 


_ B-A 
= 2 


D-C 


c = 


d = C-B, - 


B+A _D+C 

Vq ab ~ 2 ’ ^Ocd ~ 2 ’ 




c 

a 

— >0.1 
a 

0.25 

0.5 

1 

2 

OO 


i 

0.554 

0.548 

0.542 

0.535 

0.530 

0.522 

*cj 5) 

0.5 

0.542 

0.537 

0.533 

0.529 

0.525 

0.522 

K°° 

0.25 

0.532 

0.529 

0.527 

0.524 

0.523 

0.522 


0.1 

0.525 

0.524 

0.523 

0.522 

0.522 

0.522 


1 

0.554 

0.548 

0.542 

0.535 

0.530 

0.522 


0.5 

0.403 

0.393 

0.383 

0.374 

0.367 

0.361 

K* 

0.25 

0.256 

0.242 

0.231 

0.223 

0.217 

212 


0.1 

0.103 

0.091 

0.083 

0.078 

0.074 

0.072 


1 

0.600 

0.596 

0.591 

0.585 

0.581 

0.574 

K W 
K 0ab 

0.5 

0.591 

0.587 

0.583 

0.580 

0.577 

0.574 


0.25 

0.582 

0.580 

0.578 

0.576 

0.575 

0.574 


0.1 

0.576 

0.576 

0.575 

0.575 

0.575 

0.574 


1 

0.600 

0.596 

0.591 

0.585 

0.581 

0.574 


0.5 

0.470 

0.461 

0.453 

0.446 

0.440 

0.435 

K°° 

0.25 

0.343 

0.332 

0.324 

0.316 

0.312 

0.308 


0.1 

0.211 

0.202 

0.196 

0.192 

0.189 

0.187 
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Table 24: The effect of LJh on the normalized stress intensity factors at 

the maximum penetration point of two identical collinear semi- 
elliptic surface cracks under uniform tension or bending 

(Material 1) 




B+A 

Voab = ~J~ ’ 

B-A 

a ~~7 - ’ 

d = C-B, -= 1, 

C 

a 

l =1 




L 0 

— -0.2 

0.4 

0.5 

0.6 

0.8 


d 







a 







0.1 

0.832 

0.488 

0.313 

0.165 

-0.005 


0.25 

0.829 

0.481 

0.306 

0.159 

-0.006 


0.5 

0.826 

0.470 

0.299 

0.153 

-0.008 

K (#) 

K 0ab 

1 

0.821 

0.465 

0.290 

0.147 

-0.009 


2 

0.816 

0.456 

0.283 

0.141 

-0.010 


oo 

0.811 

0.446 

0.272 

0.135 

-0.011 


0.1 

0.842 

0.545 

0.397 

0.270 

0.0755 


0.25 

0.840 

0.540 

0.392 

0.266 

0.0764 


0.5 

0.837 

0.533 

0.386 

0.261 

0.0753 

K n 
0 ab 

1 

0.833 

0.526 

0.379 

0.256 

0.0742 

K°° 

2 

0.828 

0.518 

0.374 

0.252 

0.0733 


oo 

0.823 

0.511 

0.366 

0.247 

0.0725 
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Table 25: 


The effect of LJh on the normalized stress intensity factors at 
the maximum penetration point of two identical collinear semi- 
elliptic surface cracks under uniform tension or bending 

(Material 2) 


B+A 

V Oab ~ ~2~ ’ 


a = d = C-B, 



a 

k 


1 



d 

a 

T -02 

0.4 

0.6 

0.8 


0.1 

0.815 

0.444 

0.121 

-0.021 


0.25 

0.812 

0.437 

0.115 

-0.0224 


0.5 

0.809 

0.430 

0.110 

-0.0235 

K ( B ) 
K 0ab 

— 

1 

0.806 

0.423 

0.105 

-0.0245 

K°° 


2 

0.802 

0.417 

0.101 

-0.0253 


oo 

0.798 

0.409 

0.097 

-0.0260 



0.1 

0.826 

0.505 

0.234 

0.0622 


0.25 

0.823 

0.500 

0.230 

0.0611 

K (7) 
K 0ab 

0.5 

0.820 

0.494 

0.226 

0.0602 

1 

0.817 

0.488 

0.222 

0.0592 


2 

0.814 

0.483 

0.219 

0.0585 


oo 

0.810 

0.477 

0.215 

0.0578 
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Table 26 : 


The effect of LJh on the normalized stress intensity factors at 
the maximum penetration point of two identical collinear semi- 
elliptic surface cracks under uniform tension or bending 

(Material 3) 


_ B+a _ B-a 

y Oab ~ ~2 ~ ’ ° “ ~2 ~ ’ 


d=C-B, -=1, 

C 


1 



d 

a 


0.4 

0.6 

0.8 


0.1 

0.873 

0.554 

0.201 

- 0.0026 


0.25 

0.871 

0.548 

0.194 

- 0.0045 


0.5 

0.869 

0.542 

0.189 

- 0.0062 

K (- E ) 
K 0 ab 

— 

1 

0.866 

0.535 

0.183 

- 0.0078 

K°° 


2 

0.864 

0.530 

0.178 

- 0.009 


oo 

0.861 

0.522 

0.172 

- 0.0103 



0.1 

0.880 

0.600 

0.301 

0.0822 


0.25 

0.878 

0.596 

0.297 

0.081 


0.5 

0.876 

0.591 

0.292 

0.0793 

*0 aP 

— 

1 

0.874 

0.585 

0.288 

0.080 

K°° 


2 

0.872 

0.581 

0.283 

0.077 


oo 

0.869 

0.574 

0.279 

0.076 
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Table 27 


The effect of the thickness ratio on the stress intensity 
factor in a plate under uniform bending, containing a 
single crack, <r b = 6M 0 /h 2 (Material 4 and Material 5) 


a 

h 

Mat.l 

Mat.4 

Mat. 5 

0.05 

0.9885 

0.9894 

0.9878 

0.125 

0.9559 

0.9591 

0.9536 

0.25 

0.8992 

0.9059 

0.8950 

0.5 

0.8195 

0.8293 

0.8144 

1 

0.7477 

0.7571 

0.7432 

2 

0.7003 

0.7015 

0.6911 

4 

0.6701 

0.6733 

0.6675 

8 

0.6521 

0.5916 

0.5864 


*j(V 2) 
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Table 28: 


Stress intensity factors at the 
a semi-elliptic surface crack in 
tension or pure bending 


maximum penetration point of 
an infinite plate under uniform 
(Material 4 and Material 5) 


—► 0.2 


0.4 0.6 0.8 


a 

h 



0.5 

0.710 

0.303 

0.049 

-0.029 


1 

0.809 

0.441 

0.130 

-0.013 


2 

0.872 

0.575 

0.233 

0.015 

K°° 

4 

0.923 

0.704 

0.365 

0.062 


6 

0.946 

0.776 

0.458 

0.101 


Mat.4 



0.5 

0.729 

0.387 

0.173 

0.049 


1 

0.820 

0.506 

0.244 

0.071 


2 

0.881 

0.623 

0.332 

0.102 

K°° 

4 

0.928 

0.736 

0.445 

0.146 


6 

0.949 

0.799 

0.523 

0.182 


Mat. 5 


0.5 

0.697 

0.290 

0.044 

-0.0293 

1 

0.797 

0.424 

0.121 

-0.0141 

2 

0.864 

0.558 

0.221 

0.012 

4 

0.918 

0.691 

0.351 

0.057 

6 

0943 

0.766 

0.444 

0.095 


0.5 

0.717 

0.375 

0.169 

0.048 

1 

0.809 

0.491 

0.236 

0.070 

2 

0.873 

0.608 

0.322 

0.099 

4 

0.923 

0.725 

0.433 

0.142 

6 

0.945 

0.790 

0.512 

0.176 
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Figure Is 


Notation for moment and transverse shear resultants. 



2 < 



Figure 3: Geometry and notation for two collinear cracks. 
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Mat.l 

Mat.2 

Mat.3 



h/a a/h 

Figure 4: Stress intensity factor in a cracked plate under uniform bending. 
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I 

0 . 


2 . 


4 . 


a/h 

Figure 7: Stress intensity factor at the 

semi-elliptic surface crack in an 
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Mat.2 



point of a semi-elliptic surface crack in an infinite 
plate under uniform tension. 
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Appendix A 

Definition of constants A. and B j 


A 1 ~ 10k ( + 5 « S « K S 66 ) ’ 

_^ S J1 

A ’~ “i°V 

_h*_ 

A 3~ 12 / 22 ' 

h 3 2 k 

^4 = ^7 = 5 12 ) ' 

_h 2S 44 
A& ~ _10S 66 ' 

/l 2 5 55 

^6 = ~10 k 5 11 5 44 + ^12^55 " *5^ ) ’ 

/L 

j4 ® _ _ 12 / u ' 




*2 = 


*3 = 


*4 = 


*5 = 


A 3 , 

~2 A,, 

- A g , 

”^ 3^5 ’ 

“^ 4^5 _ ^ 3^6 ~ ^ 1^4 ’ 


*6 = 


^ 4^6 ^ 2^4 ^ 1^8 


*7 = 


A 2 A g . 
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Appendix B 
Expressions for R ; 

( j = 1, 2, 3 for x > 0) 

/?l( a ) = ^-( A 3 — A 2 ) [ Q 3 A 3 A 2 + ( A 3 + A 2 )a — <? 5 <* + ^6 a 1 ^i( * ) » 

* 2 ( " » = ^ A i - •'s 1 1 A A + 9 ‘ ' A . + A 3 )“ 2 - 9»* J + «X I »■( < ) ' 

« 3 ( » I = - M 1 « 3 V 2 + «,(*, + *, )« 2 - e ^ 1 + ««<>* I «,( 1 ) , 


where 

D(a) = ( Q 2 a 4 + ) |A 2 A 3 ( A 3 - A 2 ) + A,A 3 ( Aj - A 3 ) + A,A 2 ( A 2 - A, ) ] , 

m t are roots of equation (13) and A t = rrr 2 . 


Ql ^3 P l P \' 

Q 2 = p i( p s p s~W’ 

Q* = Vi . 

<?5 = V 2 - 

^6 = P 6 P 3 _ ^5^4 ' 
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P >-120 « S M S ”' 

k 3 1 

^ “ 6 S 66 

Pj = 'life V 1 S 22 S “ + S ‘ 2S,< ’ ’ 

h 5 5 55 

p < = Ims^ 3 ' 

k 3 

P 5 - 12 * S 22 > 

/» 3 1 /* 3 

k 2 k 2S < ( 

P7_ 10 « Si2S « + 10 5 66 ’ 

/i 2 /l 2 ^55 

P 8 = Ic^ 55Sl2_ los^' 
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